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We study the high-energy equivalence between helicity ±1/2 gravitinos and goldstinos in order to calculate the 
production of gravitinos in time-dependent scalar and gravitational backgrounds. We derive this equivalence for 
equations of motion, paying attention to some subtleties, mainly due to external sources, that are not present in 
the standard proofs. We also propose the Landau gauge as a simplifying alternative to the usual gauge choices, 
both for practical calculations and in the equivalence theorem proof. 



I. INTRODUCTION 



In supergravity theories |l],[2| the graviton superpartner is a spin 3/2 particle called the gravitino. This particle 
couples only with gravitational strength to the rest of matter fields, and accordingly its lifetime can be very 
long, with a decay rate of T 3 / 2 — rriy 2 /Mp. In particular, gravitinos lighter than m 3 / 2 < 100 MeV will live 
longer than the age of the Universe. This fact can have important consequences in cosmology and imposes 
stringent constraints on supergravity models. Owing to their weak couplings, gravitinos freeze out very early 
when they are still relativistic; therefore their primordial abundance can be estimated as n 3 / 2 /s — 10~ 3 
Considering only the case of unstable gravitinos, such a primordial abundance would give rise to an enormous 
amount of entropy, in conflict with the standard cosmology. In particular, gravitinos decaying during the 
nucleosynthesis can destroy the nuclei created in that era. A possible way out of this gravitino problem, is 
the existence of a period of inflation that dilutes any primordial density Q. Unfortunately the problem can 
be re-created if, after inflation, gravitinos are produced by some mechanism. In fact, this could be the case 
if during the period of inflaton oscillations, at the end of inflation, the reheating temperature was sufficiently 
high. A successful nucleosynthesis era then requires (we give some conservative bounds ||): n 3 / 2 /s < 10~ 15 for 
a gravitino mass m 3 / 2 ~ 100 GeV, n 3 / 2 /s < 10~ 14 for m 3 / 2 — 1 TeV and n 3 / 2 /s < 10~ 13 for m 3 / 2 — lOTeV. 
The production of gravitinos during reheating is due to processes involving other particles produced from the 
inflaton decay, and depends on the reheating temperature Tr as ||: n^/i/s — 10~ 14 Tr/(10 9 GeV). For a 
typical mass m 3 / 2 ~ 1 TeV, this implies Tr < 10 9 GeV. Another constraint appears in supergravity models 
where the gravitino mass is determined by the scale of supersymmetry breaking. In order to solve the hierarchy 
problem, it is then suggested that m 3 / 2 < 1 TeV pi. 

However, as a consequence of some recent works |tJ it was realized that, during the first inflaton oscillations, 
reheating can not be studied by the standard perturbative techniques. This preheating period can give rise to 
an explosive production of bosons due to the phenomenon of parametric resonance. In this period, the energy of 
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the coherent oscillations of the inflaton field is very efficiently converted into particles. In the case of fermions, 
the limit imposed by the Pauli exclusion principle avoids the explosive production, although the results still 
deviate from the perturbative expectations [|]j9|. This fact is particularly relevant when gravitinos are directly 
coupled to the inflaton, since during preheating they could be produced in excess, thus imposing new constraints 
on the particular supergravity inflationary model. 

In previous works 1 1 ll[] it was shown that the production of helicity ±3/2 gravitinos can take place during 
preheating and that the results deviate from the perturbative expectations by several orders of magnitude (see 
also fj3). In the case of helicity ±1/2 gravitinos, the production is in general more abundant, depending on 
the specific supergravity model fl3|| . Some other works dealing with this topic can be found in Jl4[ . 

In the present work we are interested in the production of helicity ±1/2 gravitinos during preheating. The 
relative difficulty of the calculations in the unitary gauge, used in the above references, suggests that we should 
explore alternative methods. In particular we will exploit the relation between helicity ±1/2 gravitinos and 



goldstinos (first pointed out in |15|) given by the Equivalence Theorem (ET) [J16|. This possibility was suggested 
in the first reference of jl3| and in the last one of |ll| it was shown how the ET could be used to study the 
helicity ±1/2 gravitino equation in the \<fi\ <C -Mp limit. 



lq ], as a way to calculate 
i are scalars and therefore 



The ET was first introduced, in the framework of non-abelian gauge theories [1^ 
processes involving longitudinal gauge bosons, but using only Goldstone bosons, whic 
much easier to handle. The first formal proof in terms of S-matrix elements was given in p3|, and that 
is basically the derivation followed within the supergravity scenario. In the last few years, and still within 
the framework of the non-abelian gauge theory, several works have completed the proof of the theorem pO , 
including renormalization effects, but also raising some questions about its Lorentz non-invariance ambiguity 
and applicability pl[ . In this paper we also discuss briefly how these new considerations may affect the gravitino- 
goldstino ET applied to the production of gravitinos during preheating. 

Intuitively, the ET tells us that, since the goldstinos disappear from the spectrum through the super-Higgs 
mechanism, giving rise to physical helicity ±1/2 gravitinos, it is possible to use goldstinos in the calculation of 
observables instead of the complicated ±1/2 gravitinos. Of course, this identification can only be carried out 
at energies high enough to neglect the masses. 

Rigorously, this theorem has only been proved for S-matrix elements containing initial or final helicity ±1/2 
gravitinos and in the absence of external backgrounds. This would provide a good approximation for gravitino 
production during the reheating period, but only at the perturbative level, where the rate of production is given 
by the decay of inflaton quanta 0,|l6| • 

However, preheating is a non-perturbative (and out of equilibrium) process and it is not obvious that the same 
proof still holds in the presence of external sources, such as the inflaton field or the space-time curvature. In 
particular, the presence of a source that creates particles makes different the initial and final vacua in the Green 
functions. In addition, these sources are present in the gauge-fixing condition, which is the starting relation in 
the ET derivations. Of course, we still expect that the intuitive relation suggested by the ET should hold, but 
since it is not the same to establish an equality at the level of matrix elements as at the operator (fields, indeed) 
level, we present in this paper a derivation more suited for the formalism in terms of equations of motion. In 
this way we can also identify the physical conditions on the sources that we need for this theorem to hold. 

Finally, we propose the Landau gauge as the best choice to perform the calculations, although, probably, it 
is not the most intuitive. In this gauge, not only the proof of the theorem, but also the final equations that 
govern gravitino production are considerably simpler. 

All the previous considerations basically concern the gravitino production process. But we also have to take 
into account the fact that we are producing very many gravitinos (out of equilibrium) which have a distribution 
in energies. Some of them will satisfy the physical conditions to apply the ET, whereas some others will not. 
Hence, we also present an additional condition on the number of those gravitinos not satisfying the applicability 
conditions, in order to obtain reliable calculations with the ET. 
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II. SUPERGRAVITY LAGRANGIAN 



Let us consider N = 1 minimal supergravity tJMJggl coupled to a single chiral superfield $, which describes 
a complex scalar field <f) and a Majorana spinor 77 satisfying tj = Cfj T = i] C , with the charge-conjugation matrix 
given by C = ij 2r y°. In principle, the derivation could be extended to more than one chiral multiplet in a 
similar way. The scalar component will play the role of the inflaton field and it will therefore be considered as 
an external background. The corresponding Lagrangian is defined by the superpotential W($) and the Kahler 
potential G($,$t) = $t$ + i og \W\ 2 . We will define: G i4> = dG/d<j>, G i4> * = dG/ d<j>* , G M . = <9 2 G/<9<W, etc. 
In this case we will have G^* = 1. The bosonic part of the Lagrangian is given by 

g- 1/2 c B = ~R + g^d^d^* + e G (3 - |G, | 2 ) , (1) 



where we are working in units M = Mp/y / 8ir = 1. In the fermionic part of the Lagrangian, we are only 
interested in those terms quadratic in the fermionic fields (gravitinos and goldstinos), since we are going to 
work with the linearized equations of motion. For the sake of simplicity we will assume that the scalar field <f> 
is real. With this assumption those terms are: 

g-^U = -l^P°^ lv D P i> a + ^fjpv + e G/2 (jfrv^r + \ {~G M - G %) m 

+ -±=G,^rr?) + -±=M^)rv, (2) 



with a^v — 5 [7^,71/]- Since we are concerned with the production of gravitinos after inflation, we assume that 
our scalar field depends only on time and that the space-time metric is of the Friedmann-Robertson- Walker 
(FRW) form. In particular, it will be very useful to work in conformal time, for which the FRW metric with 
flat spatial sections reads: 

ds 2 = a 2 {t){dt 2 -dx 2 ), (3) 



where a(t) is the Universe scale factor and the non- vanishing gravitational field is assumed to be created by the 
scalar field. 

In contrast with the ET usual proof, there are two mixing terms between gravitinos and goldstinos in eq. (Q) . 
When the scalar field has settled down at the potential minimum, <f> — ipQ, the last term does not contribute, 
and this is why it is absent from the discussions of the spontaneous breaking of supersymmetry. However, since 
we are interested in a time-dependent </>, such a term cannot be ignored any longer. In flat space-time, with 
<f> = 4>q, and when supersymmetry is not broken, i.e. G^ = 0, the mixing terms are absent and the equations 
of motion describe the gravitino evolution with only two helicity ±3/2 states. However, when supersymmetry 
is broken spontaneously, the gravitino acquires two more degrees of freedom with helicity ±1/2, because of the 
interaction with the goldstinos, giving rise to much more complicated evolution equations. 

In the unitary gauge all the goldstino dependent terms are absorbed in a redefinition of the gravitino field. 
This gauge shows explicitly the super-Higgs mechanism in which the goldstino becomes the helicity ±1/2 
components of the gravitino field. There are no mixing terms but, still, we have to deal with ±1/2 helicity 
states of a Rarita-Schwinger field, which can be rather involved. Nevertheless, the production of helicity ±1/2 
gravitinos in preheating has been calculated in the unitary gauge in pH} ]. 

Note that, although it is not necessary, we are making the inflaton responsible for supersymmetry breaking. 
Then, the inflatino also plays the role of the goldstino. This assumption simplifies the discussion since otherwise, 
and although supersymmetry would be broken during and after inflation, it would be restored at the minimum 
of the potential and the super-Higgs mechanism would not take place. Accordingly, the gravitino would not 
have a ±1/2 component. We will also assume that at the minimum the cosmological constant is zero, then 
G 2 ^ = 3G : O 0* , and we have G^ Q = \/3. 
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The equations of motion for gravitinos and goldstinos derived from eq. (g) are 

e^nD^ + l -e G ' 2 [r,l u ]^ - -J=G, e G / 2 7 ^ - -±=(^>) 7 ^ = (4) 

and 

i pr, + e G ' 2 {-G M - G%) rj - ±= e G ' 2 G^^ + -^VWW/. = (5) 

If we consider only helicity ±3/2 gravitinos, then it can be shown that the equations of motion reduce to a very 
simple form jjS: 

(ip-e G / 2 )^ 3 / 2 = (6) 

However, the helicity ±1/2 equation is much more involved and contains terms coupled to goldstinos. Nev- 
ertheless, if we are only interested in the helicity ±1/2 gravitinos high-energy behavior, E 3> w 3/ / 2 , then we 
can simplify the calculations with the ET. This limit is sensible in most supergravity inflation models with one 
chiral supermultiplct, since the typical energy of the particles created during preheating is of the order of the 
inflaton mass, m^, which is usually several orders of magnitude larger than ra 3 / 2 . For instance, in the model 
discussed in p3| , p4| ], to^ ~ 10 10 GeV, whereas m 3 / 2 < 1 TeV. Note also that all these scales are well below Mp, 
where supergravity breaks down as an effective theory. 

As it was commented before, the ET has been rigorously derived for S-matrix elements. However, to calculate 
the non-perturbative production of gravitinos during preheating, we use a formalism in terms of equations of 
motion and fields. In order to show how the gravitino-goldstino high-energy equivalence can be used in this 
context, we will follow these steps: 

i) Introduce a gauge-fixing term corresponding to a certain generalization of the gauges, which allows us 
to cancel the mixing gravitino-goldstino terms in the equations of motion. 

ii) Assume that in the asymptotic regions t — > ±oo, the external sources are static, i.e, 4> - * 4>o and <? M „ — > ?7 M „, 
and then use the equations of motion in those regions to show that d^ip^ cx m 3 / 2 ^. 

iii) Use the high-energy limit of the ±1/2 helicity projectors, -P±]y 2 = P fi / m 3/2 + £ > ( m 3/2/ E), to relate 
^±1/2 = ^±1/2^^ ^ r l' when E 3> w 3 / 2 in the asymptotic regions. 

iv) Choose the Landau gauge, £ — > oo, as an additional simplification for the calculations of goldstino 
production. 



III. GAUGE FIXING 



Goldstinos do not belong to the physical spectrum, and in the unitary gauge we can even get rid of them 
in the equations of motion. In contrast, the production of helicity ±1/2 gravitinos during reheating is gauge- 
invariant, and is only related to the goldstino production in certain gauges, called gauges, in which both 
fields appear simultaneously in the Lagrangian. Let us then consider the following gauge-fixing condition, which 



is a generalization of the R^ gauge used in [p.6 25 



When <f> is constant we recover the gauge-fixing term in Jl6[ and the limit £ — > corresponds to the unitary 
gauge. Note that in our case, due to the external sources, all the coefficients in the gauge-fixing function are 
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no longer constants. The above equation provides us with a relation between gravitinos and goldstinos, but we 
want to extract only those with helicity ±1/2, for which we will need a relation between d^ip^ an d V- I n the 
following we will use the equations of motion to obtain a relation of the desired form. 

If we assume that in the asymptotic regions t — > ±00 the space-time is flat and the scalar field settles down 
at the potential minimum (f>a, then, in those regions, the above condition reduces to 

<out w„ = yf (8) 

where m 3 / 2 = e G °/ 2 and ai„ iOIlt are the scale factor values in the asymptotic past and future. In order to simplify 
the notation, we will absorb the scale factor into the mass: m in out = a,i n ,out m 3/2] to avoid the proliferation of 
indices, we will denote m in out simply by m. With this notation, the gauge-fixing condition reads: 

wm=|v|?- (9) 



Let us recall that it is only in the static regions where the definition of particle and the separation between 
different hclicities is unambiguous. However, in the strict sense, within the inflationary cosmology neither the 
initial (t — ► —00) nor the final (t — > 00) regions can be considered static, since there is a period of inflation before 
preheating and today we know that the Universe is expanding. Nevertheless, for practical purposes, we can 
still consider the initial and final regions as static, since the particle production will mainly take place during 
the first inflaton oscillations. Accordingly, we will define our initial vacuum by imposing such initial conditions 
on our fields that they behave as plane waves before preheating. The final state has a similar behavior, since 
the rate of expansion decreases with time. Indeed, the vacuum at the end of preheating could be defined more 
rigorously as an adiabatic vacuum [ p6| , which would not yield additional gravitino production from the Universe 
expansion. 

Let us then consider first the equations of motion for gravitinos eq. and goldstinos eq. (||) in the 
initial and final regions with the notation that we have just introduced. Since the inflaton is in the minimum, 

G^*G^ 0< p = -G,0 O , and therefore 

e^ pa i 5 ^d p iPa + -m^V]^ - iJlrrvfri = 0, (10) 



/ or, — 2mi] — i\ / —rwy^t/jfj, = 0. (11) 



If we now fix the gauge using eq.(^J) in the above equations, they can be rewritten as 



e^ pCT 757^ p ^ + im[y , 7 1Vv - - 0, (12) 

3 m 2 1 

1 PV- 2mr)- 1- — -77 = 0, (13) 



In the following we will rewrite the equations of motion for goldstinos and gravitinos in the asymptotic regions 
as well as the gauge- fixing condition in a more convenient form. Contracting the gravitino equation with <9 M , 
we obtain |^7j 

im(P7"V* - l v $i>v) ~<?> (W* = 0, (14) 



whereas contracting with 7a 7^, we find 
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2;(d A 7 CT Vv- W>\) + m{nxl v ^ v + 2V-a) - 2i£ 7A ftf^ v = 0, 



(15) 



which can be contracted again with 7 A to get 

i{W>l>o - 7 A N>\) + Wl^/i - 4 ^ ft'i'^v = 0. (16) 



Substituting eq. into eq. (16), we obtain 

- i£ ^ 7 >, + \im 2 1 v i>u + 2£m = 0, (17) 

which, finally, can be rewritten as: 

(^-m + )(i^-m_)7^ = 0, (18) 

where we have defined m± — m(l ± \/l — 3/ (2£)). Note that, in the perturbative sense, the poles in the 
propagator are exactly those obtained in [ |16| . In addition we can derive the very same equation for the 
goldstino, just by multiplying eq. (|l^) by i ft: 

(i p - m+)(i p - m-)ri = (19) 

The implications of these equations are clearer if we recall that the physical fields, i.e. the ±3/2 and ±1/2 
helicity modes in the asymptotic regions, are those satisfying both 7 AI ?/v = and c^^V = 0. They have the 
correct physical mass m since they still satisfy (i ft — m)i[>^ hv = 0. In contrast, the unphysical spin-1/2 modes 
present poles in the propagator at m±, exactly as happens with the goldstinos. Therefore, by fixing the gauge 
we have only modified the poles of the unphysical modes. From ( |l4| ) and by means of the gauge-fixing condition, 
we obtain 

^m(2 frfipv - 2d v ip u ) - i ft^rrnj = (20) 

Since r/ satisfies eq. (Jis|) , we have two possible solutions: (i ft — m+)r) — and (i ft — mJ)r\ = 0; together with 
the gauge-fixing condition (p0|), they yield 



r,- 9^-^-771^ = 0. (21) 



From this expression we get: 



^-l/lf('-^) 



(22) 



At first sight, this equation relates the unphysical gravitino d^ip^ with the goldstino; however, the key obser- 
vation is that, as we will show, at high energy, d^ipn tends to the physical helicity ±1/2 gravitino. Note that, 
apparently, there are two relations, one for goldstinos that correspond to the m_ solution and another for those 
with m_|_. 



IV. THE EQUIVALENCE THEOREM 



In the asymptotic initial and final regions, we expect that a general solution of the equations of motion for 
gravitinos and goldstinos will be written as a linear superposition of on-shell positive and negative frequency 
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plane waves |27f| . In particular, let us consider a negative frequency mode solution in the initial region with 
momentum p^ — [u>i n ,0, 0,p) (there is no loss of generality in this choice), with p^p^ = mf n and p = \p\, such 
that p 3> rriin : 

= 3/2 ' ^Mp) + O (^f) , (23) 

where V'm(p) ^ s the corresponding Fourier component. For the unphysical goldstino we have 

V P (x) = 3/a ] e^iiip) + O (^) . (24) 

Similar expressions can be written for the positive-frequency solutions and for solutions in the final region. Note 
that the on-shell conditions for gravitinos and goldstinos are different because of the different positions of the 
poles. In particular, for physical gravitinos we have p^p^ — mf n and for goldstinos and unphysical gravitinos 
(7 M VV an d d^tpfj,) we have p^ — (m™) 2 (we will not use m_ for reasons that will become clear below). We 
have thus included the 0{m in /p) term at the end of (|2|) and ([m|). Strictly, it should not be there for physical 
gravitinos, but at this level we keep a compact notation between physical and unphysical gravitinos. 

The spin-1 polarization vectors are given by e fl (p, m) — a(t)5®e a (j), m), where 

e a (p,l) = A=(0,l,i,0) ,e a (p,0) = —(p,0,0,Uin) ,e„(?,-l) = -- ^ (0,1, -*,0). (25) 

V 2 ITT-in V 2 

If u(p,s) are spinors with definite helicity s — ±1/2, then P±u(p, ±1/2) = u(p, ±1/2), where P± = (1/2)(1 ± 
l5l^^tj,(p, 0)) are the helicity projectors. Accordingly, the helicity ±3/2 and ±1/2 projectors are nothing but 

P± 3/2 = P±e M (p,±l), 

If 1/2 = \[\ P ^M ±1) + \J\P±*M 0)- (26) 



We see that, at high energy , the ±1/2 projector behaves as 

3 m m \ p 



P± 1/2 = \^P±^ + o(^), (27) 



where we have neglected e M (p, ±1) with respect to e M (p, 0). Let us then define the helicity ±1/2 components of 
the gravitino field in the asymptotic initial regions and in momentum space as 



^±1/2 (?) = P±i/ 2 ^n(P) = 



6 mi. 



O 



m i7 
V 



Vv(£)- 



(28) 



At high energies, we see that the helicity ±1/2 gravitino tends to the unphysical d^ipu field and therefore we 
can use the gauge-fixing condition in ( p2| ) to obtain a relation between each Fourier mode of the goldstino and 
the helicity ±1/2 gravitino. (As pointed out in p6| , it is essential that both d^tp^ and rj have the same poles, 
in order to rewrite the equality in eq. ( |22] ) in terms of Fourier modes) . 

For arbitrary values of the £ parameter it would be necessary to take into account both solutions, with m+ 
and with m_. In this case, the solution of the goldstino equation is 

r) P {x)=r]t{x)+Vp{x). (29) 
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From (P4), we see that each solution is related to the 1/2 helicity gravitino with different proportionality 
constants, i.e. for negative frequency solution we will have in the initial region 



^±1/2 (?) = 



i-e- 



p. 



±1/2 



o 



v 



fj + ' (p). 



(30) 



In the S-matrix derivations of the ET, the proportionality constant between the helicity ±1/2 gravitinos and 
the goldstinos disappears once the external lines of the Green functions have been removed, the momenta are 
on-shell, and the tensor indices are contracted with the corresponding polarization vectors [l6| . However, this is 
not so straightforward in the "semiclassical" proofs based on the generating functional formalism, either within 
supergravity jl6| or even in the non-abelian context |l7j. These "semiclassical" proofs are given for the clever 
choice £ = 3/2, where to_ = m + and the proportiona ity constant is unique (in the non-abelian case the choice 
is £ = 1 and the proportionality constant is unity). 

However, for our purposes, it is much more appropriate to choose the Landau gauge. Indeed, in an arbitrary 
generalized i?j gauge, eq. (^) will be written 



c.G/2 



(G ^ + G%) r, - e G / 2 G^^e G ' 2 G.^ = 



(31) 



The presence of the last term makes it very difficult to obtain solutions even numerically. However, we get 
a dramatic simplification by using the Landau gauge, £ — > 00, in which the last term, which is the most 
complicated, vanishes. Thus we have: 

i p v - e G ' 2 {G M + G%) V = (32) 

Note that this last expression corresponds to the m + case in (fl9|) . The m_ = solution is just an artifact due 
to the multiplication by i $ and we do not have to consider it in the previous formulae. Therefore, eq.(p8|) now 
reads 

v (f). (33) 



4>±1/2(P) 



2iP ± + — 

p 



This is the relation we were looking for. Note that this is an equality at the level of fields and not for S-matrix 
elements. This result shows that although the helicity ±1/2 gravitinos and the goldstinos can evolve differently 
during the oscillations of the scalar field, they approach each other in the asymptotic regions (up to a constant). 
The fact that this result is valid only in the asymptotic regions, is sufficient for our purposes. Since ([52|) is 
just the standard equation of motion for a fermion field in curved space-time with a time-dependent mass, it is 
straightforward to apply the standard techniques of particle production. 

Once more we stress that such a result is only useful when the energy of the particles we are producing is 
much larger than their masses. In our case, the inflaton is in a frame where it is homogeneous, only depends 
on i, and oscillates with a typical frequency m$. Thus, we expect the physical momentum of the gravitinos to 
be 0(m<p), so that we can use eq. (|33| ) if m 3 / 2 <C m^,. Usually we will evaluate a out right after the preheating 
ends and then mi n ^ out <C m 3 / 2 - The fact that we are in an homogeneous background is technically relevant due 
to the remarks about the ET and Lorentz invariance done in pl| in the context of non-abelian gauge theories. 
Indeed the ET is not Lorentz invariant, since not only is the helicity decomposition frame-dependent, but also 
terms like 0(m/E), 0(m/p)... do have very different values depending on the reference frame. The fact that 
the inflaton field is homogeneous ensures that all our gravitinos are produced from rather similar conditions 
and we can apply the ET to the vast majority of them. 
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V. PARTICLE PRODUCTION 



Up to now our discussion has been purely classical. In order to interpret these solutions in terms of particle 
number, we have to quantize them p6| , |28| ], which has already been done in the unitary gauge |l3| . We are 
interested in the Landau gauge, £ — > oo, where the goldstino equation of motion (|32] ) reduces to that of a 
Majorana fermion coupled to a scalar field in a curved space-time, whose quantization is also a well known 
problem of particle production J§|,|9| . 

Let us then consider a classical solution to eq.(|32|) with helicity I, such that in the past (t — > — oo), it behaves 
as a negative-energy plane- wave, i.e.: 

rjf(x) - \ e^-Vutf, I), (34) 

a in \J 2uJ tn 

where a in is the scale factor at the end of inflation. In the asymptotic future (f — > oo), because of the presence 
of the time-dependent background fields, this solution will no longer behave as a negative-energy mode; rather, 
it will be a linear superposition of positive and negative frequency modes 

r$(x) - 3 1 (a; ^^Mp, + P1 P ^-^-^(-p, I)) , (35) 

where, for a given p, we have {wf n out ) 2 — (mf n out ) 2 +p 2 . The Bogolyubov coefficients satisfy 

K/ + |$/ = l- (36) 

Using again our previous result in (|3^) we can identify each Fourier mode above with the corresponding Fourier 
mode for a helicity ±1/2 gravitino (up to a constant). In particular we will find that the Bogolyubov coefficient 
for the helicity ±1/2 gravitino will be the same as that of the goldstino, up to 0{m/p), i.e f?r l = (l+O(m/p))0^ t . 
Notice that because of the different masses of goldstinos and physical gravitinos, the correction to the Bogolyubov 
coefficients for gravitinos can depend on time as exp(iA^i), where 0Jg Ut = u> ou t + A w . However, such a term 
will be relevant only when t ~ 1/A. U , that is, much later than the end of preheating. Moreover, remember 
that the Bogolyubov coefficients are normalized according to ( |36| ) and therefore the proportionality constant is 
irrelevant. 

As a remark, let us note that if we had a renormalizable theory whose low-energy limit is supergravity, we 
could still use our estimates with the ET, irrespective of the renormalization corrections [ p0[ needed in the 
complete proof of the theorem. The reason is that we do not need the proportionality constant to obtain the 
Bogolyubov coefficients. 

Therefore the number of gravitinos created with helicity I = ±1/2 and momentum p, will be given by 

Ki? (37) 

In conclusion, solving the equation of motion for the goldstinos in the presence of the external backgrounds, 
and using (|37|), we will obtain the number of helicity ±1/2 gravitinos created during preheating. In the above 
expression it is explicit that only the knowledge of the solutions in the asymptotic regions is relevant to the 
particle number calculation. 

As an analytic check of the ET, we can compare with the unitary gauge results of |l3| , obtained in the global 
supersymmetric limit, \<j)\ <C Mp, We only have to note that, when \<p\ <C Mp, our goldstino equation in the 
Landau gauge ([52]) is reduced to 
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i fl V - {d^W) V = 0, 



(38) 



which is the very same equation obtained in [|f3|] for the helicity ±1/2 gravitinos. Note however, that the 
condition \(f>\ <C Mp, is not necessary to prove the ET. Our high-energy result is valid independently of the size 
of |0|. 



VI. NUMERICAL EXAMPLE 



As a possible application of the previous results and for further comparison with other works |l3j, we will 
study a simple supergravity model based on the superpotential 

W = VX— . (39) 



Unfortunately, at the minimum of the corresponding inflaton potential, supersymmetry is restored. We will 
then assume that ( |39"| ) is only valid far from the minimum and that close to = it is modified to satisfy 
the assumptions needed to apply the ET (see paragraph before eq.(|J)). For this superpotential, the effective 
gravitino mass oscillations are not damped in time. Still, we can apply the ET if we take cf)(t) = for t < 
and t > nT, where n is an integer (in our example n = 6) and T is the inflaton oscillation period. Despite 
these problems, we consider this model useful as a numerical illustration of the ET. Hence, we have taken the 
initial amplitude of the inflaton oscillation to be <po — Q.2Mp (where </> = \/2Fte0 is the canonically normalized 
inflaton), which implies that the effective goldstino mass is oscillating with amplitude ma — V2XM. With 
these initial and final conditions we have calculated numerically the number of goldstinos produced from eq. 
( pj^ ) , using the standard results for the production of fermions obtained in || . We thus look for solutions of the 
goldstino equation with momentum p and helicity I of the form 



rjP'ix) = a- 3/2 (tyP- s UP l (t), 



(40) 



with 



U*(t) = 



1 



*7% - V ■ 7+ a(t) (e G/2 {G M + G 2 ) )] f p i(t)u(p, s) 



(41) 



Using the above ansatz, we can write the equation as follows: 
— + k 2 + — ^(he G ' 2 (G ^ 



G 



—e G (G.d>d> + G 2 j.) 



A 



f K i(t) = 



(42) 



with k = p/ (ainVX) and t = ai n \f\t and the new scale factor is defined as b(t) — a(t)/a,i n . The initial conditions 
are f K i(0) = 1 and f K i(0) — —in. In particular, for the goldstino occupation number we find 

N v Kl (nT) = i- (2k + i[f* Kl (nT)f Kl (nT) - f* Kl (nT) } UnT)]) (43) 

Using eq.(|37|), we get the occupation number of helicity ±1/2 gravitinos directly from ([i"3|), whereas for those 
with helicity ±3/2 we use eq. (||) following the same steps as before jio]]. In Fig. 1, we have plotted both 
spectra and we see that for this particular model, the production of helicity ±3/2 gravitinos is suppressed by 
two to three orders of magnitude with respect to the ±1/2 gravitinos. Note that Nki depends on the number of 
oscillations and on the inflaton initial conditions. Thus, we cannot make a comparison for the whole spectrum 
with the unitary gauge estimations obtained in the first references in [O and the complete numerical analysis 
of the last reference in jl3| (made without the ET, although they have checked that the formulas hold in one 
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particular case). Nevertheless, the ET tells us that, since the mass of the goldstino will be less than a(t)>/X, both 
results should agree for p » a(nT)\/X, i.e. with the above definitions k 3> b(nT) (in our example b(nT) ~ 12), 
irrespectively of the initial conditions. Indeed, the orders of magnitude are in good agreement with the previous 
works. 




However, there are some other models in which the application of the ET is somewhat restricted, such as a pure 
quadratic supcrpotcntial W = m^$ 2 . In this case, at the minimum of the potential, supersymmetry is restored 
and therefore the definition of helicity ±1/2 gravitinos is meaningless. Ignoring this problem, although the 
amplitude of the inflaton oscillations is damped, the goldstino mass term, which contains the second derivative 
of the superpotential, tends to a constant of 0{m^), typically much larger than m 3/ / 2 - Then, the ET will be 
useful to calculate only that portion of the spectrum with energy much higher than the inflaton mass. 

Finally, there is an additional condition to obtain reliable predictions. Since gravitinos are produced in large 
numbers, with vastly different energies, we have to estimate what fraction of them does not satisfy the ET 
conditions. For that purpose we define the number density of gravitinos produced with both helicities as 

1 f°° 

nit)= ^mL N - ip2dp - ^ 

Then, for those gravitinos with momenta lower than their mass, i.e. p 2 < a 2 ut m 3/2' ^ nc ET does not apply. 
Thus the number density of "excluded" gravitinos is 

n{t)exc < 3l 2a3( t ) tt ^ m 3/2 ( 45 ) 

In addition, we can estimate the total number density of produced gravitinos as: 

" (<)tot ~ 3 ^a3 W Q m w ( 46 ) 

Hence we obtain that the fraction of gravitinos that do not satisfy the ET applicability conditions is 

n e xc < / aoutm 3/2 \ 3 
ntot ~ V a in m.^ / 

Accordingly, the additional ET applicability condition is: a utm 3 / 2 « a^m^. We see that the result depends 
on the duration of the preheating era and the ratio of gravitino and inflaton masses. Typically, the production 
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takes place in a few inflaton oscillations, which implies that the scale factor only grows by a few orders of 
magnitude, not enough to overcome the mass difference. Therefore, in these models, the ET safely describes 
the production of the vast majority of helicity ±1/2 gravitinos. 



VII. CONCLUSIONS 



We have studied the production of helicity ±1/2 gravitinos using the equivalence of goldstinos and gravitinos 
at high energies. We have shown that in the gauges, where the goldstino and gravitino equations of motion 
decouple, the classical solutions of the goldstino equation are proportional to those of the helicity ±1/2 gravitinos 
in the asymptotic static regions. This result is sufficient to relate the production of goldstinos to the production 
of helicity ±1/2 gravitinos. Furthermore, we have shown how in the Landau gauge, the equation of motion for 
the goldstino is considerably simpler. As a check, we have compared our results with previous ones obtained in 
the unitary gauge and we have found good agreement in the equivalence theorem applicability regions. 

We have clearly identified the ET applicability conditions in the context of gravitino production: i) The 
frequency of the inflaton field oscillations should be larger than the gravitino mass, 3> m 3 / 2 - If we are 
interested in the pure gravitational production, then one should also require H 3> "1.3/2- These conditions 
ensure that the typical energy of the produced particles will be larger than their masses ii) The sources should 
vanish asymptotically, which implies that the space-time curvature should decrease with time and also that the 
amplitude of the inflaton oscillations should be damped, iii) The ET calculations are only useful if most of the 
gravitinos have a large enough energy, which requires a out m 3 / 2 "C a^m^. 

Concerning the potential applications of our results, we have shown that cubic superpotentials, with a slight 
modification in their form, together with appropriate initial conditions, could satisfy the above applicability 
requirements. More elaborated models would require to extend the present proof to more than one chiral 
multiplet and non-minimal supergravity. Further work along this line is in progress. 
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